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The Potts model is one of the most popular spin models of statistical physics. Prevailing majority of the work 
done so far corresponds to the lattice version of the model. However, many natural or man-made systems are 
much better described by the topology of a network. We consider the g-state Potts model on an uncorrelated 
scale-free network for which the node-degree distribution manifests a power-law decay governed by the ex- 
ponent A. We work within the mean-field approximation, since for systems on random uncorrelated scale-free 
networks this method is known often to give asymptotically exact results. Depending on particular values of q 
and A one observes either a first-order or a second-order phase transition or the system is ordered at any finite 
temperature. In a case study, we consider the limit q = 1 (percolation) and find a correspondence between 
the magnetic exponents and those describing percolation on a scale-free network. Interestingly, logarithmic 
corrections to scaling appear at A = 4 in this case. 
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A paper dedicated to Mykhajlo Kozlovskii 
on the occasion of his 60th birthday. 

1. Introduction 

Considerable attention has been paid recently to analysis of phase transition peculiarities on 
complex networks [l|, Q ■ Possible applications of spin models on complex networks can be found in 
various segments of physics, starting from problems of sociophysics [3j to physics of nanosystems [J], 
where the structure is often much better described not by geometry of a lattice but by a network. In 
turn, the Potts model, being of interest also for purely academic reasons, has numerous realizations, 
see e.g. (jjj for some of them. The Hamiltonian of the Potts model that we are going to consider in 
this paper reads: 

- H = 2 X! J iJ S ^i^j + X! h rfm,o, (1-1) 

here, q is the number of Potts states (ni = 0, 1, ...q — 1), hi is a local external magnetic field chosen 
to favour the 0-th component of the Potts spin variable n^. The main difference with respect to the 
usual lattice Potts Hamiltonian is that the summation in (jl.ll) is performed over all pairs i, j of TV 
nodes of a network, being proportional to the elements of an adjacency matrix of the network. 
For a given network, equals J if nodes i and j are linked and it equals otherwise. 

Being one of possible generalizations of the Ising model, the Potts model possesses a richer 
phase diagram. In particular, either first or second order phase transition occurs depending on 
specific values of q and d for d-dimensional lattice systems Q . It is well established by now that 
this picture is violated by introducing structural disorder, see e.g. Q and references therein for 2d 
lattices and for 3d lattices. Here, we will analyze the impact of changes in the topology of the 
underlying structure on thermodynamics of this model, when Potts spins reside on the nodes of 
uncorrelated scale-free network, as explained in more details below. 
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The notion of complex network which is actively exploited in the physical community already 
for several decades, from the mathematical point of view is nothing else but a random graph 
Accordingly, analysis of the Potts model on various graphs has already a certain history. Similar to 
the Ising model [&] , the Potts model on a Cayley tree exhibits no longe-range order [9( but exhibits 
a phase transition of continuous order [l(| • Some exact results for the three state Potts model with 



competing interactions on the Bethe lattice are given in and the phase diagram of the three 



states Potts model with next nearest neighbour interactions on the Bethe lattice is discussed in 



12J, Potts model on the Apollonian network was considered in |13| . So far, not too much is known 



about critical properties of the model (jl.ll) on scale free networks of different types. Two pioneering 
papers [3 [l| (the latter paper was further elaborated in [l6j]) that used the generalized mean 
field approach and recurrent relations in the tree-like approximation, respectively, although agree in 
principle about suppression of the first order phase transition in this model for the fat-tailed node 
degree distribution, but differ in the description of the phase diagram. Several MC simulations also 
show evidence of changes in the behavior of the Potts model on a scale-free network in comparison 



with its 2d counterpart ly, |17(. The Potts model on an inhomogeneous annealed network was 
considered in relation between the Potts model with topology-dependent interaction and 
biased percolation on scale- free networks was considered in |l_9] . Rounding of the first order phase 
transition for the Potts model with large values of q on scale-free evolving networks was observed 
in [H. 

In this paper, we will calculate thermodynamic functions of the Potts model on an uncorrelated 
scale-free network. In contrast to [3], we will work with the free energy (Ref. [l4| dealt with the 
equation of state). This will enable us to get a comprehensive list of scaling exponents governing 
second-order phase transition as well as percolation exponents. We will show the emergence of log- 
arithmic corrections to scaling for percolation and calculate the logarithmic correction exponents. 
The paper is organized as follows. In section [5] we derive general expressions for the free energy 
of the Potts model on uncorrelated scale-free network. Thermodynamic functions will be further 
analyzed in section [3l where we obtain leading scaling exponents and show the onset of logarithmic 
corrections to scaling for some special cases. In section [4] we will further elaborate the q = 1 limit 
of the Potts model, that corresponds to percolation on a complex network. We end by conclusions 
and outlook in section [5] 

It is our pleasure to contribute by this paper to the Festschrift dedicated to Mykhajlo Kozlovskii 
on the occasion of his 60th birthday and doing so to wish him many more years of fruitful scientific 
activity. 

2. Free energy of the Potts model on uncorrelated scale-free network 

In the following we will use the mean field approach to analyze thermodynamics of the Potts 
model (fTTTjl on an uncorrelated scale-free network, that is, a network that is maximally random 
under the constraint of a power-law node degree distribution: 

P(k)=c x k~ x , (2.1) 

where P(k) is the probability to find a node with degree k and c\ can be readily found from the nor- 
malization condition y~], —u P(k) = 1, with k+ and k* being the minimal and maximal node degree, 
correspondingly. For an infinite network, limyv-i-oo k* — >• oo. A model of uncorrelated network with 
a given node-degree (called also configuration model, see e.g. [2l[) provide a natural generalization 
of the classical Erdos-Renyi random graph and is an undirected graph maximally random under 
the constraint that its degree distribution is a given one. It has been shown, that for such networks 
the mean field approach leads in many cases to asymptotically exact results. In particular, this has 
been verified for the Ising model using recurrence relations 1221 and replica method [26| and further 
applied to 0(m)-symmetric and anisotropic cubic models [23[, mutually interacting Ising models 
24 1 as well as to percolation For the Potts model however, two approximation schemes, the 



mean field treatment [lj] and an effective medium Bethe lattice approach [15[ were shown to lead 
to different results. 
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2.1. General relations 

To define an order parameter and to carry out the mean field approximation in the Hamiltonian 
(jl.ip . let us introduce local thermodynamic averages: 



Mi — fintfi , v i — Sni,ayiQ: (2-2) 

where the averaging means: 

Sp(...)exp(-ff/T) 



(•••)= ^••• / J V ' \ (2-3) 

T is the temperature and we choose units such that the Boltzmann constant fcg = 1. The partition 
function 

Z = Spexp(- J ff/T), (2.4) 



and a trace is defined by: 



N q-l 

SP(...)=I]^(-)- (2-5) 



i=l n;=Q 



The two quantities defined in (|2.2I) can be related using the normalization condition S nu o + 
Eo=i S ni,a = I, leading to: 

Vi = {l-m)/{q-l). (2.6) 

Observing behaviour of averages (|2.2[) calculated with the Hamiltonian in the low- and high- 
temperature limit: fj>i(T — > oo) = Ui(T — > oo) = l/q, ^i(T — >■ 0) = 1, Vi(T — > 0) = the local order 
parameter (local magnetization), < m, < 1, can be written as: 

mi = — . (2.7) 



Now, neglecting the second-order contributions from the fluctuations S ni , nj — ^n;,n.,- one gets the 
Hamiltonian in the mean field approximation: 

- H mfa = JijSntfimj + - Y Jiji 1 ~ m 3 + (! - timmj) + ^ ^5 ni ,o- (2.8) 

The free energy in the mean field approximation, —g = TlnSpcxp(— H m f a /T), readily follows: 

1 ■ J- m • + h' 
-9=~J2 J «( 1 - m i + C 1 - + T E ln t ex P ( — 'V' ! ) + « - !]• ( 2 - 9 ) 

As usual within the mean field scheme, the free energy (12.91) depends, besides the temperature, 
both on magnetic field and magnetization. The latter dependence is eliminated by the free energy 
minimization, leading in its turn to the equation of state. For the Potts model on uncorrelated scale- 
free network the equation of state, that follows from (|2.8[) was analyzed in Ref. [14]. Here, we aim 
to further analyze temperature and field dependency of the thermodynamic functions. Opposite 
to the mean field approximation for lattice models, where one assumes homogeneity of the local 
order parameter (putting m., = m for lattices) , intrinsic heterogeneity of a network, where different 
nodes may have in principle very different degrees, does not allow to make such an assumption. One 
can rather assume within the mean field approximation that the nodes with the same degree are 
characterized by the same magnetization. Therefore, the global order parameter for spin models 
on network is introduced via weighted local order parameters (see e.g. 2JJ). Following [14| let us 
define the global order parameter by: 

m = (2.10) 
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Within the mean field approach, we substitute the matrix elements Jij in (|2.9I) by the probability 
Pij of nodes i, j to be connected. The last for the uncorrelated network depends only on the node 
degrees fcj, kf 

J« = JPii = Jjfifc, (2-11) 

where J is an interaction constant, (k) — l/N^2iLi is the mean number of linksQ The free 
energy (12.9[) . being expressed in terms of (|2.10l) . (12.11[) . contains sums of unary- type functions 
over all network nodes. Using the node degree distribution function, these sums can be written as 

sums over node degrees: jj Yli=i f(^i) = J2k=k P{k)f(k)- In the infinite network limit, N — > oo, 
k* — > oo, passing from sums to integrals and assuming homogeneous external magnetic field hi = h 
we get for the free energy of the Potts model on an uncorrelated scale-free network: 

nn 

Jk + Jk m+ Jk(q-l) m2 _ rln(e ^$±A + g _ !)] P[k)dk f (2.12) 

K 

where the node-degree distribution function is given by (12.11) . For small external magnetic field h, 
keeping in (I2.12[) the lowest order contributions in h, hm and absorbing the m-independent terms 
into the free energy shift we obtain: 

oo 



g = — ^-m + '-m z - T / ln(e mJfc/J +q- l)P(k)dk - ^ ' rah . (2.13) 

q q J q 2 T 

Free energy (|2.13[) is the central expression to be further analyzed. In spirit of the Landau theory, 
expanding (|2.13l) at small m and first keeping terms ~ m 2 one gets for the above expression at 
zero external magnetic field: 

- Inq + (T - T )m 2 , (2.14) 

where To = 2q(fc) • Provided that the second moment (fc 2 ) of the distribution (|2.ip exists, one 
observes that depending on temperature T, the coefficient at m 2 changes its sign at To. This 
temperature will be further related to the transition temperature. Another observation, usual 
for the spin models on scale free networks Q is, that the system remains ordered at any finite 
temperature when (k 2 ) diverges (since To — > oo). For distribution (|2.1j) this happens at A < 3. 
Therefore, we will be primarily interested in temperature and magnetic field behaviour of the 
Potts model at A > 3|j The expansion of the function under the logarithm in (|2.13p at small 
order parameter m involves both the small and the large values of its argument, mJk. To further 
analyze (|2.13p , let us rewrite it singling out the contribution (|2.14[ 3 and introducing new integration 
variable x — mJk/T: 

J(k)(q-l) 2 c^J)" J (g-i)(fc)J ( 

(T-To)m + — / (p(x)dx — mh , (2.15) 



3 qT T x 

where a;* = mJk*/T and 



<p(x) = [-ln(e x +q-l)+lnq+- + %-^-x 2 ]\. (2.16) 



2<j 2 



1 Such approximation makes the model alike the Hopfield model used in description of spin glasses and autoas- 
sociative memory |2Sl430ll . 

2 Scale- free networks with fc* = 1 do not possess a spanning cluster for A > A c (A c = 4 for continuous node degree 
distribution and A c ~ 3.48 for the discrete one [3lll. We can avoid this restriction by a proper choice of fc* > 1. 
3 Again, we absorb the constant — lnij into the free energy shift. 
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Note that the Taylor expansion of the expression in square brackets in f|2 . 16[) at small x starts from 
a; 3 , whereas at large x the function ip(x) behaves as x 2 ~ x and therefore the integral in (|2.15p is 
bounded at the upper integration limit at A > 3. To analyze the behaviour of the integral at the 
lower integration limit when m-^Owe proceed as follows. 



2.2. Non-integer A 

Let us first consider a case when A is non-integer. Then we represent ip(x) for small x asQ 



[A-l] 



(2.17) 



i=3 



=w 



where [£] means an integer part of i, a, = ai (q) are the coefficients of the Taylor expansion: 

oo 

- ln(e x + q - 1) = ^ a i x i , (2.18) 

The first coefficients are as follows: 



i=0 



a = -lng, a\ = —1/q, a 2 = 



g-1 

2q 2 : 



a 3 = 



(g-l)(g-2) _ 
6q 3 



04 = 



(q-l)(q 2 -6g + 6) 
24g 4 



This is an integration of the first sum in (12.171) that leads to terms that diverge at x — > 0. Let us 
extract them from the integrand and evaluate the integral in (I2.15[) as follows: 



lim / cp(x)dx = lim / <p(x) — > 
:*^0 / z«— >o / L ^ — ' 



[A-l] 



[A-l] 



dx 



i=3 



lim > 



i=3 



-dx . 



(2.19) 



For the reasons explained above, the first term in (|2.19p does not diverge at small to anymore, 
neither it diverges at large x, so one can evaluate this integral at to = numerically. We will 
further denote it as: 

[A-l] 



c(q,X) 



i=3 



dx . 



(2.20) 



Numerical values of c(q,X) at different q and A are given in Table [TJ 
Integration of the second term in (|2.19p leads to: 



[A-l] J a [A-l] 
i=3 „ i—3 



CLi(x+) 



-A+i+1 



A— 1 — « 



(2.21) 



Finally, substituting (|2.20p and (|2.2ip into (|2.15p we arrive at the following expression for the first 
leading terms of the free energy at non-integer A: 



J(k}(q-1) 2 , c x c(g,X) A _j aijmJk+y x _j 

" = — (T-T )m + yA „ 2 (toJ) + ca 2^ A _ _ T - 

i=3 



J(fc)(«-i; 



q 2 T 



mh + 0(m lx] ) 



(2.22) 



4 It is meant in 112. 17H and afterwards, that the first sum is equal to zero if the upper summation limit is smaller 
than the lower one, i.e. for A < 4. 
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Table 1. Normalized numerical values of the coefficient c(q,X)/(q — 1), Eq. ()2.20|) . for different 
q and A. 



X\q 


1 


2 


3 


4 


6 


8 


5.4 


-3.0692 


-0.0079 


-0.0002 


0.0013 


0.0011 


0.0007 


5.1 


-5.9318 


-0.0454 


-0.0106 


-0.0006 


0.0028 


0.0025 


4.8 


-0.1686 


0.0352 


0.0148 


0.0058 


0.0005 


-0.0005 


4.5 


-0.2439 


0.0237 


0.0154 


0.0085 


0.0030 


0.0012 


4.2 


-0.6809 


0.0275 


0.0344 


0.0240 


0.0119 


0.0067 


3.9 


0.5975 


0.0420 


-0.0540 


-0.0528 


-0.0346 


-0.0231 


3.6 


0.7240 


0.0830 


0.0065 


-0.0076 


-0.0102 


-0.0085 


3.3 


1.4001 


0.2469 


0.0790 


0.0315 


0.0052 


-0.0010 



2.3. Integer A 

Let us consider now the case of integer A. To single out the logarithmic singularity in the 
integral of Eq. (I2.15[) let us proceed as follows [32j . Denoting 



K(y) = j tp(x)dx (2.23) 
y 

we take the derivative with respect to y: 

d -^ = -<p(y). (2.24) 
ay 

Now, K(y) can be obtained expanding the expression in square brackets in (|2.16[) at small y and 
integrating Eq. (HOI)) : 

/ oo — A 

Viv)dy= £ T-i-l - a *-iMv)+C(q,\), (2.25) 

i=3,i/A-l 

with an integration constant C(q, A) and coefficients dj given by (|2.18[) . Numerical values of C(q, A) 
at different q and A are given in Tabic [2] 

Substituting K(mJk+), cf. Eq. (I2.23[) . into (|2.15p we arrive at the following expression for the 
free energy at integer A: 

J(k)(q - 1) 2 CaQa-1 , n A-li , \ru \\ 

g = — (T-To)m - rA _ 2 (mJ) mm + c\[C{q, A) — 

«A-t ln(^/T)]i^ + c A g ^^T^ - J -^^mh + 0(ro W) (2.26) 



i=3 



Expressions (|2.22[) . (|2.26l) for the free energy of the Potts model will be analyzed in the subse- 
quent sections in different regions of q and A. 



Table 2. C(q, X)/(q - 1) for different q and A. 



A\g 


1 


2 


3 


4 


6 


8 


4 


0.9810 


0.0355 


-0.0085 


-0.0134 


-0.0109 


-0.0079 


5 


0.4853 


0.0194 


-0.0527 


-0.0483 


-0.0303 


-0.0197 
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3. Thermodynamic functions 

Having in mind analysis of the Potts model also in the percolation limit q = 1, let us rescale 
the free energy by the factor (q — I): g m f a = g/(q — 1) and absorb it by re-defining the free energy 
scale. Then, each term in (I2.22p . (|2.26[) is also to be divided by (q — 1). Let us use the following 
notations for several first coefficients at different powers of m in (|2.22l) . (j2.26[) : 



.4 
B 
C 
K 
D 



2J{k) 

q 

CA (JfcQ 3 (g-2) 

2<7 3 (A-4) ' 
c x (Jk*) 4 (q 2 -6g + 6) 
6q 4 (A - 5) 



B 1 



c X J 3 {q-2) 



c 



2q 3 

c A J 4 (g 2 -6g + 6) 
6q 4 



cxJ^ciq.X) 



(9-1) 



J(k) 



(3.1) 
(3.2) 
(3.3) 
(3.4) 
(3.5) 



Below, we will start analysis of the thermodynamic properties of the Potts model by determining 
its phase diagram in different regions of q and A. 



3.1. The phase diagram 

To analyze the phase diagram, let us write down expressions of the free energy at small values 
of m, keeping in (|2.22p . (I2.26[) only the contributions that on the one hand allow to describe non- 
trivial behaviour, and, on the other hand, ensure thermodynamic stability. Since the coefficients 
at different powers of m are functions of q and A, cf. (|3.1[) - (|3.5I) . the form of the free energy will 
differ for different q and A as well. 




T<Tc T<Tc 

T=Tc T=Tc 

T>Tc T>Tc 

a b 



Figure 1 . Typical behaviour of the free energy of the Potts model on uncorrelated scale free net- 
work at zero external field h = 0. a: continuous phase transition; b: first-order phase transition. 
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3.1.1. 1 < q < 2 

As far as the coefficients K, B and B' at ra A_1 , m 3 and m 3 lnm are positive in this region of 
q, it is enough to consider only the three first terms in the free energy expansion: 

3<A<4: g = A( T - T )m 2 + ^m^ 1 - ^mh, (3.6) 

A = 4: g = ^-{T - T Q )m 2 + -^m 3 In - - ^mh, (3.7) 
21 61 z m 1 

A>4: g = A.( T -T )m 2 + ^m 3 -^mh. (3 . 8) 

Typical m-dependence of functions (|3.6[) ~ (13.8[) at h = is shown in Fig.[T^. As it is common for the 
continuous phase transition scenario, the free energy has a single minimum (at m = 0) for T > Tq. 
A non-zero value of m that minimizes the free energy appears starting from T = Tq. In particular, 
the transition remains continuous in the percolation limit q = 1, as will be further considered in 
sections l3~!Q HI 

3.1.2. q = 2 

For q = 2, the Potts model corresponds to the Ising model. Indeed, in this case the coefficient 
at m 3 vanishes and the first terms in the free energy expansion read: 

3<A<5: g= A( T - T Q )m 2 + K^rn^ 1 - ^mh, (3.9) 
A = 5: g= A( T _ To )m 2 + ^m 4 In 1 - ^m/i, (3.10) 
A>5: 5 = A(T_ To)m 2 + _^ TO 4_^ m/l . (3 . n) 

It is easy to check that the above coefficients if, C, C" are positive for q — 2. Therefore, again the 
free energy behaviour corresponds to the continuous second-order phase transition, see Fig. [1^. 

3.1.3. q > 2 

In this region of q, phase transition scenario depends on the sign of the next-leading contribution 
to the free energy. Indeed, for positive K the free energy reads: 

3<A<A C ( (Z ): g = A. {T ^ To )m 2 +K 1 ^m x - 1 -^mh, (3.12) 

where K remains positive in the region of A bounded by the marginal value A c defined by the 
condition 

c{q,\ c ) = 0, (3.13) 

with c(q,X) given by (|2.20p . The free energy p,12j) is schematically shown in Fig. [T£i for different 
T. As in the former cases, 13.1.11 13.1.21 it corresponds to the continuous phase transition. With an 
increase of A, the coefficient K becomes negative and one has to include the next term: 

A C ( 9 )<A<4: 5 = A( T _ To ) TO 2 + ^_L_ m A-i + _|_ m 3_^ m ^ ( 3 . 14) 

B > for A < 4. Now, because of the negative sign of the coefficient at m A_1 the free energy 
possesses a local minimum at a decrease of T (see Fig. [TJa) and the order parameter manifests 
a jump at the transition point T c : scenario, typical for the first order phase transition. With 
further increase of A, one has to include next terms in the free energy expansion for the sake 
of thermodynamic stability. However, the sign at the second lowest order term remains negative, 
which corresponds to the free energy behaviour shown in Fig. [TJd: the phase transition remains of 
the first order. 
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The above considerations can be summarized in the "phase diagram" of the Potts model on 
uncorrelated scale-free network, that is shown in Fig. O There, we show the type of the phase 
transition for different values of parameters A and q. 




Figure 2. The phase diagram of the Potts model. The black solid line separates the 1 st order 
phase transition region from the 2 nd order phase transition region (shaded). The critical ex- 
ponents along the line are A-dependent. In the 2 nd order phase transition region, the critical 
exponents are either A-dependent (below the light solid line, red online) or attain the mean field 
percolation values (above the red line). For q — 2, A > 5 (shown by the black dotted line), the 
critical exponents attain the mean field Ising values. Two different families of the logarithmic 
corrections to scaling appear: at A = 5, q — 2 (a square, red online) and atA = 4, l<g<2 
(light solid line, red online). For A < 3 (the region below the dashed black line) the system re- 
mains ordered at any finite temperature. The values for the rest of critical exponents are listed 
in Tables El 11 



3.1.4. Vq, 2 < A < 3 

As it was outlined above, for 2 < A < 3 the Potts model remains ordered at any finite temper- 
ature. Similar as for the Ising model 22], it is easy to find the high-temperature decay of the order 
parameter in this region of A for any value of q > 1. Since (fc 2 ) becomes divergent for 2 < A < 3 
one does not single out this term in the expression for the free energy (cf. (|2.f 5p ). As a result, the 
corresponding expressions for the free energy read: 

2 < A < 3 : g = ^m 2 + ^m x ~ l - ^mh, (3.15) 

» o A 2 C A J 2 2 1 C A J 2 r C'((Z,3) I Jk* , 2 D 

A = 3 : q = —=m H In =— » ln( — — m — —mh . 

y 2T Aq 2 T m T 1 q - 1 2q 2 T J T 

(3.16) 
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Here, the expressions for the coefficients K' , C'(q, 3) are the same as for K, C(q,3), Eqs. Q3.4p . 
(|2.25[) with the only difference that the function f(x) used for their calculation does not contain 
the x 2 term. It is easy to check that the free energies p,15[) . p. 161) are minimized at any finite 
temperature by the non-zero value of m that decays at high T as [14l . |26| : 



2 < A < 3 
A = 3 



rn 
m 



Te~ aT , 



a > 0. 



(3.17) 
(3.18) 



The above equations (|3 . 1 7(1 , (|3.18[) give the temperature behaviour of the mean- field order param- 
eter m. The connection with the magnetization M is found from the self-consistency relation: 



M 



(dg_ 
\dh 



(3.19) 



One can check, that the solution of this equation at large T is of the form M ~ ^ . Correspondingly, 
this leads to the following high temperature decay of M 22|: 



2 < A < 3 
A = 3 



M 
M 



e~ aT 1 a>0. 



(3.20) 
(3.21) 



For the sake of simplicity, in what follows below we will express the thermodynamic functions 
in terms of the mean field order parameter m. To get their M dependence, one has to take into 
account the above considerations. 



3.2. Regime of the second order phase transition, critical exponents 

Let us find the critical exponents, governing behavior of thermodynamic functions in the vicinity 
of the 2nd order phase transition point r = T — Xo|/Xb = 0, h = 0. To this end, we will be interested 
in the following exponents, that govern temperature and field behavior of the order parameter m, 
the isothermal susceptibility xt — (^)t, the specific heat Ch = T(§^)h, and the magnetocaloric 
coefficient itit = — T"(ff )t- 

h = 0; m~T?, Xt^t-"*, Ch ~ T - a , m T ~ t~ u , (3.22) 
t = 0: to-Zi 1 / 5 , xt~/i~ 7c , c h ~br a % m T ~h- u '. (3.23) 

As in the former subsection, we analyze this behaviour in different regions of q and A. The results 
of this analysis are summarized in Table [3l 



Table 3. Leading critical exponents of the Potts model on uncorrelated scale-free network. 



q 


A 


a 


a c 


P 


5 


7 


7c 


to 


UJ C 


1 < q < 2 


3 < A < 4 


A-5 

A-3 


A-5 

A-2 


1 

A-3 


A-2 


1 


A-3 
A-2 


A-4 

A-3 


A-4 

A-2 


1 < q < 2 


A > 4 


-1 


-1/2 


1 


2 


1 


1/2 








q = 2 


3 < A < 5 


A-5 
A-3 


A-5 

A-2 


1 

A-3 


A-2 


1 


A-3 
A-2 


A-4 
A-3 


A-4 
A-2 


q = 2 


A > 5 








1/2 


3 


1 


2/3 


1/2 


1/3 


q>2 


3 < A < X c (q) 


A-5 
A-3 


A-5 
A-2 


1 

A-3 


A-2 


1 


A-3 
A-2 


A-4 
A-3 


A-4 
A-2 



3.2.1. 1 < q < 2 

In this region of q, the free energy is given by the expressions (|3.6[l - (|3.8[l . For T > To, g is 
minimized at h = by a zero value of the order parameter m = 0. For T < T$, the minimum 
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of the free energy corresponds to the non-zero m. Based on the expressions (|3.6I) - (|3.8I) we find in 
different regions of A: 

3<A<4: m= r ^(_£_)^, (3 . 24) 
AT 2 

A = 4: m = -^Lr\\nr\-\ (3.25) 
AT 2 

A>4: m= —^t. (3.26) 
B 

By formulas (13.241) . (13.261) at q = 1 we reproduce the corresponding results for the percolation on 
scale-free networks [33j: the usual mean field percolation result for the exponent j3 = 1 for A > 4 
and /J = for 3 < A < 4. Note the appearance of the logarithmic correction at the marginal 
value A = 4. The obtained values of the exponent are given in Table [3J Subsequently, we obtain 
the rest of the exponents defined in (|3.22[) . (|3.23[) and display them in the first two rows of Table 
[3] as well. 

Similar to the order parameter, (|3.25p . the temperature and field behaviour of the rest of 
thermodynamic functions at A = 4 in the vicinity of the critical point is characterized by the 
logarithmic corrections. Let us define the corresponding logarithmic-correction-to-scaling exponents 
by 0: 

h = 0: m ~ tP\ lnr f, X T ~ t~t| lnr c h ~ t~ q | lnr f , m T ~ t^I lnr |", (3.27) 

t = 0: m ~ h 1 ' 5 ] In ft \\ X t ~ ^~ 7c | In ft | fc , c h ~ ft- Qc | In h f c , m T - fr u ' \ In h f° . 

(3.28) 

For their values we obtain the numbers given in Table[U Note, that all exponents are negative: loga- 
rithmic corrections speed up the turning to zero of the decaying quantities and weaken singularities 
of the diverging quantities. We discuss this behaviour in more details in Section 2J 



Table 4. Logarithmic-corrections exponents for the Potts model on uncorellated scale-free net- 
work. 



q 


A 


a 


d c 


P 


5 


1 


7c 


Co 




l<q<2 


A = 4 


-2 


-3/2 


-I 


-1/2 





-1/2 


-1 


-1 


q = 2 


A = 5 


-1 


-1 


-1/2 


-1/3 





-1/3 


-1/2 


-2/3 



3.2.2. q — 2, the Ising model 



For different A, the free energy is given by (|3.9[) - (|3.11[) . Minimizing these expressions one finds 
for the order parameter at h = 0, T < T : 



3 < A < 5 



-L n 



\K(X-1)J 



m 



ATI 



C" 



•Vaiinrl-i/a, 



A > 5 



AT a r i/2 
C 



(3.29) 
(3.30) 
(3.31) 



Corresponding critical exponents for the other thermodynamic quantities are given in the third 
and fourth rows of Table [3] Logarithmic corrections to scaling (13.27[1 - f|3 . 28[) appear at A = 5, their 
values are given in the second row of Table |H Critical behaviour of this model on an uncorrelated 
scale-free network was a subject of intensive analysis, see e.g. the papers 22|, |23|, [25j, l26( and by 
the result given in Table [3] we reproduce the results for the exponents obtained there. 
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3.2.3. q > 2, 3 < A < A c (g) 

In this region of q, X the phase transition remains continuous, see the phase diagram, Fig. [21 and 
the free energy is given by the expression p. 121) . Correspondingly, one finds that the spontaneous 
magnetization behaves as 

3<A<A c (g): m = Tp{ ^-^ )^'. (3.32) 

The values of the rest of the critical exponents are given in the sixth row of Table [3] Since the 
leading terms of the free energy (|3. 121) at 3 < A < X c (q) coincide with that of the Ising model at 
3 < A < 5, (|3.9[) . the behaviour of thermodynamic functions in the vicinity of the second order 
phase transition is governed by the same set of the critical exponents: the Potts model at q > 2, 
3 < A < X c (q) belongs to the universality class of the Ising model at 3 < A < 5. This result was 
first obseved in [l4| by treating the mean field approximation for the equation of state. 



3.3. The first order phase transition 

For q > 2, A > X c (q) the phase transition is of the first order, see the phase diagram on Fig. 
O As we have outlined in Sec. 13.1. 3[ the next-leading order term of the free energy has a negative 
sign and the free energy behaves as shown in Fig. [TJj. As further analysis shows, the higher the 
value of A the more terms one has to take into account in the free energy expansion in order to 
ensure the correct g(m) asymptotics. Therefore, in the results given below we restrict ourselves 
by the region X c (q) < A < 4, where the free energy is given by Eq. (|3.14p . The first order phase 
transition temperature T c is found from the condition <?(m = 0,T C ) = g(m ^ 0,T C ), see the red 
(middle) curve in Fig.[T)3: 

r.-r„ + f .„. -* J * "<*-W-«> , (3.33) 



.3K(X-1)(X-3)J A 
For the jump of the order parameter Am at T c we find: 

Am=( , ~ B , -)^~\ (3.34) 

\3K(X-1)(X-3)J y ' 

Another thermodynamics function to characterize the first order phase transition is the latent heat 
Q. It is defined by: 

Q = AS ■ T c , (3.35) 

where AS is the jump of entropy at T c . With the free energy given by (I3.14[) we find the entropy 
as 

Considering the entropy at the transition temperature we can find the latent heat at the first order 
phase transition for X c (q) < A < 4: 

Q = 4(Am) 2 . (3.37) 



4. Notes about percolation on scale-free networks 

By the results of section 13.2.11 we cover also the case q — 1, that corresponds to percolation 
on uncorrelated scale-free networks. The 'magnetic' exponents governing corresponding second 
order phase transition are given in Tables G2 HI Let us discuss them in more details, in particular 
relating them to percolation exponents. The following exponents are usually introduced to describe 
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behavior of different observables near the percolatiorJl point p c (40l | : the probability that a given 
site belongs to the spanning cluster 

Poo ~ (P-Pcf, P>Pc, (4.1) 

the number of clusters of size s 

(4.2) 



-r -si a* 



S e 



\P-Pc\-", (4.3) 



the cluster size at criticality 
the average size of finite clusters 

(s)~\p- Pc \~\ (4.4) 

The above defined exponents /3 and 7 coincide with the 'magnetic' exponents j3 and 7 of the q = 1 
Potts model (see Tables [21 SJ) ■ Therefore, the probability that a given site belongs to the spanning 
cluster and the average size of finite clusters for percolation on uncorrelated scale-free networks 
are governed by the scaling exponents: 

Ht 3 a>:, 4, <«> 

7 = 1, A > 3. (4.6) 

The exponents r and a may be derived with the help of familiar scaling relations (40| : 

ap = t - 2, (4.7) 

7 = (3 - t)/<t. (4.8) 

Substituting the values of /3 and 7 (|4.5p . (|4.6p into (|4.7|) . (|4.8[) one arrives at the following expres- 
sions for the exponents r and a: 



2A-3 
A-2 ' 
5 

2 - 

A-3 
A-2 ' 



( , 9 , 

From the high-temperature behaviour of the Potts model magnetization at 2 < A < 3, (|3.20[) . 
one arrives to the scaling exponents j3 = 1/(3 — A), 7 = —1 for the corresponding observables for 
percolation at p c = 0. 

By formulas (|4.5[) . (14.61) . (|4.9p . and (|4.10l) we reproduce the results for the scaling exponents 
that govern percolation on uncorrelated scale- free networks 3^, 34 1 and also those that were found 



for the related models of virus spreading [35|, [36J. All the above mentioned papers do not dis- 
cuss explicitly the case A = 4 and possible logarithmic corrections that arise there. Moreover, a 
recent review ^37], that also discusses the peculiarities of percolation on uncorrelated scale-free 
networks does not report about logarithmic corrections (its Eq. (95) is perhaps wrong since it 
gives no logarithmic corrections for A = 4.) By our results in the first row of the Table 2] we give 
a comprehensive list of critical exponents that govern logarithmic corrections to scaling appearing 
for the Potts model at q = 1, A = 4, as correctly predicted within the general Landau theory 



for systems of arbitrary symmetry on uncorrelated scale-free networks |41| . One may compare 
our values with the corresponding exponents of the d-dimensional lattice percolation at d = 6: 
a = j3 = 7 = 5 = a c — 2/7 (see e.g. |38j). In this respect the logarithmic-correction exponents 
for the lattice percolation at d = 6 and for the scale-free network percolation at A = 4 belong to 
different universality classes. It is easy to check, that quoted in the last row of Table S] exponents 
obey the scaling relations for the logarithmic-corrections exponents: $(5 — 1) = 88 — 7, a = 2/3 — 7 

1-2) (4- 
/3+7 



5 For definitcncss, let us consider the site percolation and denote by p here and below the site occupation proba- 
bility. 
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5. Conclusions and outlook 

In this paper we have analyzed the critical behaviour of the g-state Potts model on an un- 
correlated scale-free network. The mean held approach we use in our calculations often leads to 
asymptotically exact results when critical behaviour on an uncorrelated scale-free network is con- 
sidered. However, in the case of the Potts model, two similar approximate schemes of calculations, 
the mean held [lj| and the recurrent relations for the tree-like random graphs [l5[ differ in their re- 
sults concerning the phase diagram. In particular, for q > 2 and A < 3 both approaches predict that 
the system always remains ordered at finite temperature, however for A > 3 depending of specific 
values of A the first approach predicts the first or the second order phase transition, whereas the 
second approach predicts the first order phase transition. Our results complete the above analysis, 



however, unlike the Ref.[l4|, where the equation of state was considered, we have considered the 
thermodynamic potential which enabled us to present comprehensive analysis of temperature and 
magnetic field dependence of thermodynamic quantities. It is worth to note here, that the model 
considered here is alike the Hopfield model [28[ , for which the mean field approximation is known 
to give exact results. 

Our main results are summarized in Fig. [5] and in Tables 131111 Depending of the values of q and 
of the node degree distribution exponent A, the Potts model manifests either the first-order or the 
second-order phase transition or it is ordered at any finite temperature, see Fig. [51 In the second 
order phase transition region (shaded in the figure), it belongs either to the universality class of 
the Ising model on an uncorrelated scale-free network (with A-dependent critical exponents) or it 
is governed by the mean field percolation (1 < q < 2, A > 4) or mean field Ising (q — 2, A > 5) 
exponents. 

One of the major differences of the critical behaviour of the Potts model in the second order 
phase transition regime in comparison to the Ising model that its logarithmic corrections exponents 
belong to two different universality classes. As it is well known, in certain situations the scaling 
behaviour is modified by multiplicative logarithmic corrections (see (38| for a recent review and 



42j for the specific case of the Potts model). For lattice systems, such corrections are known to 
appear, in particular, at the so-called upper critical dimension, above which the mean field regime 
holds. For the scale-free networks, where the very notion of dimensionality is ill-defined, a change 
in the node degree distribution exponent A may turn a system to such regime. For the Ising model 
it happens at A = 5 and is caused by the divergencies in the fourth moment of the node degree 



distribution (k A ) 1 as it was analyzed in details in |23|. In addition to these corrections, for the 
Potts model we observe an onset of multiplicative corrections to scaling at 1 < q < 2 for A = 4. 
Contrary to the Ising case, these are caused by the divergencies in the third moment of the node 
degree distribution (fc 3 ). This difference in the reason of their appearance causes also the difference 
in their numerical values: this new set of the logarithmic correction to scaling exponents belongs 
to the new universality class. In particular, they govern percolation on uncorrelated scale free 
networks at A = 4. 
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KpHTHHHa noBe^iHKa Mo^ejii IIoTTca Ha cKjia^HHx Mepe»cax 



iHCTHTyT 4>13HKH KOHfleHCOBaHHX CHCTGM HAH yKpaiHH, UA- 79011 JIbBIB, YKpaiHa 

iHCTHTyT >K. Jiainypa, YHiBepcHTeT JIoTapHHfii', F-54506 BaH^yBp jie HaHci, ©paHiria 

Mo^ejiB noTTca OflHa 3 HaimonyjisipHiHiHx MO^ejieii CTaTHCTHHHoi <pi3HKH. BijiBinicTB po6iT, BiiKOHaHiix 
paHime, CTOcyBajiacB rpaTKOBOi Bepcii uiei MO^ejii. OflHaK, 6araTO npnpoflHHx Ta CTBopeHHx jho^hhoio 
CHCTeM Ha6araTO Kpame onncyiOTBca Tonojiorieio MepejKi. Mh po3rji5iflaeMO g-CTaHOBy MO^ejiB noTTca Ha 
HecKopejiBOBaHift 6e3MacniTa6HiH Mepe>Ki i3 CTeneHeBO 3racHOio (pyHKuieio po3no^rijiy CTyneHiB By3jiiB i3 
noKa3HiiKOM A. npainoeMO b Ha6jin>K6HHi cepe^HBoro nojiH, ockIjibkh ^jth chctcm Ha HecKopejiBOBaHHx 
6e3MacniTa6HHx MepejKax pen MeTOfl nacTO ,a,03BOji5ie otphmsth acHMnTOTHiHO TOHHi pe3yjiBTaTH. B 
3ajie>KHOCTi Bi^ 3H&MGHL q Ta A cnocTepiraeMO cpa30Bi nepexo^H nepnioro ^,pyroro po,zry, a6o jk 
CHCTeiua 3ajiHHiaeTBCH BnopaflKOBaHOio npn 6yn,B-HKifl TeiunepaTypi. TarcojK po3rjiHflaeMO rpaHHino 
q = 1 (nepKOjiHuia) Ta 3naxo^HMO Bi^noBiflHicTB Mi>K MarHiTHiiMH kphthmhhmh noKa3HHKaMH Ta 
noKa3HHKaiviH, mo onncyiOTB nepKonsnriio Ha 6e3MacTa6mH Mepejid. HiicaBO, mo b HBOMy BHnaflKy 
jiorapncpMiHHi nonpaBKH flo CKeftjimry 3'hbjihiotbch npn A = 4. 

KjHOHOBi cjiOBa: Mo^ejiB noTTca, CKjiaflHi MepejKi, nepKOjismia, KpHTH<mi noKa3HHKH 




PACS: 64.60.ah, 64.60.aq, 64.60.Bd. 
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